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Sergio Colafrancesco, Stefano Profumo, '0 and Piero Ullio 3 >0 

1 INAF - Osservatorio Astronomico di Roma, via Frascati 33, 1-00040 Monteporzio, Italy and 
Istituto Nazionale di Fisica Nucleate, Sezione di Roma 2, F00133 Roma, Italy 
2 Division of Physics, Mathematics and Astronomy, California Institute of Technology, 
Mail Code 106-38, Pasadena, CA 91125, USA 
3 Scuola Internazionale Superiore di Studi Avanzati, Via Beirut 2-4, 1-34014 Trieste, Italy and 
Istituto Nazionale di Fisica Nucleare, Sezione di Trieste, 1-34014 Trieste, Italy 

We explore the possible signatures of dark matter pair annihilations in the nearby dwarf spheroidal 
galaxy Draco. After investigating the mass models for Draco in the light of available observational 
data, we carefully model the dark matter density profile, taking advantage of numerical simulations 
of hierarchical structure formation. We then analyze the gamma-ray and electron/positron yield 
expected for weakly interacting dark matter particle (WIMP) models, including an accurate treat- 
ment of the propagation of the charged particle species. We show that unlike in larger dark matter 
structures - such as galaxy clusters - spatial diffusion plays here an important role. While Draco 
would appear as a point-like gamma-ray source, synchrotron emission from electrons and positrons 
produced by WIMP annihilations features a spatially extended structure. Depending upon the cos- 
mic ray propagation setup and the size of the magnetic fields, the search for a diffuse radio emission 
from Draco can be a more sensitive indirect dark matter search probe than gamma rays. Finally, 
we show that available data are consistent with the presence of a black hole at the center of Draco: 
if this is indeed the case, very significant enhancements of the rates for gamma rays and other 
emissions related to dark matter annihilations are expected. 

PACS numbers: 95.35.+d,12.60. Jv,98.70.Rz,98.56.Wm 
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I. INTRODUCTION 



The astrophysical search for signals of dark matter (DM) particle pair annihilations in cosmic structures on large 
scales (from galaxies to clusters of galaxies) is, potentially, a very powerful technique, highly complementary to direct 
DM searches, in the quest for the identification of the fundamental nature of DM. The widest and more definite set of 
results can be harvested through a multi- frequency survey of DM annihilation signals over the whole electromagnetic 
(e.m.) spectrum (see, e.g. [jj, hereafter CPU2006, and references therein) by using a detailed treatment of both 
the microscopic interaction properties of the hadronic and leptonic secondary yields of WIMP annihilations, and of 
the subsequent emissions originating by the yields themselves in the astrophysical environment at hand. Various 
astrophysical systems have been taken into consideration to this aim. The central regions of ordinary galaxies (like 
our own Galaxy) are usually considered among the best places to set constraints on the presence and on the nature 
of DM particles (see, e.g, Q for a review, and the analyses in [J|j][I|- among others). 

The typical faintness of DM signals within viable WIMP scenarios makes, in fact, the Galactic Center, or the 
central regions of nearby galaxies (like M31), the most plausible and promsing places to detect signals of WIMP 
annihilations. However, the expected DM signals have to contend, there, with the rich and often poorly understood 
astrophysical context of thermal and non-thermal sources (SN remnants, pulsars, molecular clouds, to mention a few), 
whose spectral energy distributions (SEDs) cover the whole e.m. spectrum, reaching even TeV energy scales, as the 
recent results from HESS, MAGIC, Cangaroo have clearly shown (see, e.g., |(J and references therein; see, however, 
also 7]). In this respect, galaxy cores are likely not the best places to definitely identify DM annihilation signals. 

Galaxy clusters have the advantage to be mass-dominated by DM and, in some cases, like the nearby Coma cluster, 
to have a quite extended spectral and spatial coverage of thermal and non-thermal emission features which enable to 
set interesting constraints on the properties of DM (see, e.g., 0,001 and refs. therein for various aspects of the DM 
SEDs in clusters). The study of the DM- induced SEDs in galaxy clusters has been shown to be quite constraining for 
DM WIMP models, and can even be advocated to shed light on some emission features (e.g., radio halos, hard- X-ray 
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and UV excesses, gamma-ray emission) which are still unclear. Nonetheless, the sensitivity and spatial resolution of 
the present and planned experiments in the gamma-rays, X-rays and radio do not likely allow to probe more than a 
few nearby clusters. It is therefore, mandatory to remain within the local environment to have reasonable expectations 
to detect sizable emission features of possible DM signals. 

Globular clusters have also been proposed (see e.g., [10() as possible sources of gamma rays from WIMP annihila- 
tions, but with expected signals well below the sensitivity threshold of future experiments, mainly due to their quite 
low mass-to-light ratios. 

The ideal astrophysical systems to be used as probes of the nature of DM should be mostly dark (i.e., dominated 
by DM), as close as possible (in order to produce reasonably high fluxes), and featuring central regions mostly devoid 
of sources of diffuse radiation at radio, X-rays and gamma-rays frequencies, where the DM SEDs peak (see, e.g. 
CPU2006 for general examples). 

Dwarf spheroidal (dSph) galaxies closely respond to most of these requirements, as they generally consist of a stellar 
population, with no hot or warm gas, no cosmic-ray population and little or no dust (see, e.g., for a review). 
Several dSph galaxies populate the region around the Milky Way and M31, and some of them seem to be dynamically 
stable and featuring high concentrations of DM. 

Among these systems, the Draco dSph is one of the most interesting cases. This ob j ect has already been considered 
as a possible gamma-ray source fed by DM annihilations in recent studies [l^, Ha. fl4L fl5L Hq | , in part triggered by an 
anomalous excess of photon counts from Draco reported by the CACTUS collaboration in a drift-scan mode survey 
of the reg ion surrounding the dSph galaxy [l7| . The nature of the effect is still controversial, but it has been shown 
in [T3.[l5| to be in conflict, in most WIMP models, with the EGRET null-result in the search for a gamma-ray source 
from the direction of Draco Ha|. Other gamma-ray upper limits have been obtained by the Whipple 10-m telescope 
collaboration as well (see e.g., [19|V 

The observational state-of-the-art for Draco goes, however, beyond gamma-ray emissions: radio continuum upper 
limits on Draco have been obtained by Fomalont et al. [2(j with the VLA. These authors report an upper limit of 
J v < 2 mJy at v = 4.9 GHz (this is a 3a level limit). Typical magnetic field strength of B ~ 2 — 4 fiG for dwarf 
galaxies similar to Draco have also been derived from radio observations at 5 GHz [2l| . The X-ray emission from the 
central part of Draco has an upper limit provided by ROSAT [22j. The count rate detected by the PSPC instrument 
in the (0.1-2.4) keV energy band is < 0.9 • 10~ 3 s~ 1 corresponding to an unabsorbed flux limit of Fx < 1.7 • 10 -14 
erg cm _2 s _1 . This flux corresponds to an X-ray luminosity upper limit of Lx < 0.01 • 10 36 erg s^ 1 . 

The main point we wish to make in the present analysis is that a complete multi- frequency analysis of the astrophys- 
ical DM signals coming from Draco might carry much more information, and can be significantly more constraining, 
in terms of limits on DM WIMP models than, for instance, a study of the emissions in the gamma-ray frequency 
range alone. 

As we show in the present analysis, available observational data, and the possible detection of WIMP annihilation 
signals from Draco by future instruments can be, in principle, of crucial relevance for the study of the nature of 
WIMP DM: the expected emission features associated to DM annihilation secondary products are, in fact, the only 
radiation mechanisms which can be expected in a system like a dSph, as originally envisioned by Colafrancesco [9l.l23|. 
Following our original suggestions, and pursuing the systematic approach we outlined for the case of Coma (see CPU 
2006), we present here a detailed analysis and specific predictions for the WIMP DM annihilation signals expected 
from Draco in a multi- wavelength strategy. 

Specifically, we first derive the DM density profile of Draco in a self-consistent ACDM scenario in Sect. II. We 
then discuss the gamma-ray emission produced in Draco from DM annihilation, assuming a set of model-independent 
WIMP setups |lj. Gamma-ray emissions, and constraints, are studied in Sect. III. We then present in Sect. IV the 
signals expected from Draco at all frequencies covered by the radiation originating from the secondary products: 
synchrotron emission in the radio range, Inverse Compton scattering of electrons and positrons produced by DM 
annihilation off CMB and starlight photons, and the associated SZ effect. We also discuss, in Sect. V the possible 
amplification of these signals by an intervening black hole at the center of Draco. We present our conclusions in the 
final Sect. VI. 

Throughout this paper, we refer to the concordance cosmological model suggested by WMAP 3yr. 24] ; namely, 
we assume that the present matter energy density is ft m = 0.266, that the Hubble constant in units of 100 km s _1 
Mpc -1 is h = 0.71, that the present mean energy density in baryons is Of, = 0.0233//i 2 , with the only other significant 
extra matter term in cold dark matter He dm — H m — ^f>, that our Universe has a flat geometry and a cosmological 
constant A, i.e. = 1 — fl m , and, finally, that the primordial power spectrum is scale invariant and is normalized 
to the value as = 0.772. 
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II. THE DARK MATTER DENSITY PROFILE IN DRACO 



Modeling the distribution of dark matter for dSph's is not a straightforward task. The radial maps of the star 
velocity dispersions clearly indicate that dSph are dark matter dominated systems. However, available observational 
data do not provide enough information to unequivocally determine the shape and concentration of the supporting 
dark matter density profiles (see e.g. the recent analysis of Ref. [25| for the case of Draco, under investigation here). 
Such freedom is partially reduced restricting to ACDM inspired scenarios, as appropriate for dark matter in the form of 
cold WIMP particles. Within this structure formation picture, numerical N-body simulations of hierarchical clustering 
predict that Milky Way size galaxies contain an extended population of substructures, with masses extending down 
to the free streaming scale for the CDM component (as small as 10~ 12 — 10~ 3 Af Q in the case of neutralinos in 
supersymmetric models or in other WIMP setups |2fiLl27j). and surviving, at least in part, to tidal disruption: dwarf 
satellites stand as peculiar objects, since they are the smallest ones featuring a stellar counterpart, while mechanisms 
preventin g st ar formation are supposed to intervene for lighter objects (among scenarios supporting this interpretation, 
see, e.g., [28]). In case of isolated CDM halos, properties of the dark matter density profile have been investigated 
in some detail through numerical simulations: a universal shape and a correlation (on average) between the object 
mass and its concentration are expected (more details will be given in the following section) . The picture is less clear 
for satellites, like Draco, standing well within the dark matter potential well of the hosting halo. Tidal forces may 
have significantly remodeled the internal structure of these objects, an effect which is likely to depend, e.g., upon the 
merging history of each satellite. Based again on numerical simulations, significant departures from the correlation 
between mass and concentration parameter observed for isolated halos have been reported in the literature, as well 
as discrepant results regarding whether the universal shape of the density profile is preserved [2^| or not|3(| in the 
subhalos, after tides have acted and these systems have reached a new equilibrium configuration. 



A. Mass models within the ACDM framework 



The main dynamical constraint we consider for mass models for the Draco dSph is the observed line-of-sight 
velocity dispersion of its stellar population. The underlying, necessary, assumption we shall make here is that the 
stellar component is in equilibrium, and hence that the Jeans equation applies to this system; if this is the case, one 
finds that the projection along the line of sight (l.o.s.) of the radial velocity dispersion of stars can be expressed in 
terms of M(r), the total (i.e. including all components) mass within the radius r (' [31II32I ): 



where v(r) is the density profile of the stellar population and S(i?) represents its surface density at the projected 
radius R. In the derivation of Eq. Q), we have assumed that the anisotropy parameter is constant over radius; 
in terms of the radial and tangential velocity dispersion, respectively a r and erg, = 1 — a^/a^: = 1 denotes the 
case of purely radial orbits, = that of a system with isotropic velocity dispersion, while — > —00 labels circular 
orbits. As we will see shortly, the anisotropy parameter is important since we recover in our analysis the well known 
degeneracy between the reconstructed mass profile and the assumed degree of stellar anisotropy. 

Extensive photometric studies are available for the Draco dwarf; we refer to the analysis in Odenkirchen et al. [3i| 
relying on multicolor data from the SDSS (sample and foreground determination labeled S2 in that analysis) and 
reproduce the result for the radial profile of the surface brightness in Fig. ^ (left panel) . We also show two alternative 
fits of the data: one option is the generalized exponential profile proposed by Sersic [34| and implemented in the case 
of Draco also by Lokas, Mamon & Prada [35| : 

E(#) = £ exp[-(i?/iis) 1/m ], (2) 

choosing the parameter 1/m = 1.2, and fitting the scale radius Rs and central surface brightness So to the data (the 
best fit procedure gives Rs — 7'3) . As second possibility, we follow Mashchenko et al. [23 and consider a modified 
Plummer model: 

E(fl) = E„[l+(J2/iip) 2 ]~ (o ~ 1)/2 , (3) 

setting the exponent a = 7, and then fitting the value for the scale radius (Rp = 14f6) . For each of the two E(i£), 
the luminosity density profile v(r) is obtained by inverting the definition of surface brightness with the Abel integral 
formula, i.e. implementing the de-projection: 
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FIG. 1: Left panel: radial profile for the surface brightness distribution of stars in Draco; data are from Odenkirchen et al. |33| 
while fits are with a Sersic profile or a modified Plummer model. Right panel: the corresponding luminosity density profiles. 



The inversion is performed numerically for the Sersic profile, while it can be done analytically for the modified Plummer 
model; results are shown in Fig. ^ (right panel) and one can see that the mild differences in the surface brightness 
are only marginally amplified in the luminosity density profiles. Here we are referrin g to luminosities in the V-band 
and, following a gain [35| . we have adopted for the distance of Draco the value 80 kpc [2f|, or, equivalently, a distance 
modulus of 19.5 |32|, standing in between (and in agreement at 1 a) the other recent estimate of 75.8 ± 0.7 ± 5.4 kpc 
from Ref. [3^| and the value of 82 ± 6 kpc from the compilation of Mateo . To add the stellar component in the 
total mass term M(r) in Eq. QJ, we need an estimate for the stellar mass-to-light ratio in the V-band; we mainly 
refer to one of the largest values quoted in the literature, Ty = 3M q /Lq, including in it a possible subdominant gas 
component. 

The ansatz we implement for the dark matter component is that of a spherical distribution sketched by a radial 
density profile: 



P(r) = p'g(r/a) ; 



(5) 



given in terms of the function g(x) and of two parameters, i.e. a scale radius a and a normalization factor p' . This is 
in analogy with the usual description of dark matter halos from results of numerical N-body simulations in terms of 
a universal density profile; we take as guideline for our mass models the form originally proposed by Navarro, Frcnk 
& White E3: 



9nfw(x) 



1 



x(l + x) 2 



(6) 



and a shape slightly more singular towards the center proposed by Diemand et al. [41| (hereafter labeled as D05 
profile): 



gD05(x) 



1 



As a further option, we consider the Burkert profile (|4^): 

9b(x) 



with 7 ~ 1.2 . 



1 



(l + x)(l + x 2 ) ' 



(7) 



(8) 



i.e. a model with a large core, in agreement with the gentle rise in the inner part of the rotation curves occurring in a 
vast class of galaxies, including dwarfs Mechanisms of gravitational heating of the dark matter by baryonic 
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FIG. 2: Line of sight velocity dispersion for the three best-fit-models considered in the paper. Data are from Munoz et al. 2005 
(black filled circles) and from Wilkinson et al. (red filled squares). 



components during or after the baryon infall have been advocated to reconcile these observations with the central 
density cusps of the profiles introduced above [51 14fll4n1 | : these models are still contrived and it is probably premature 
to say whether in the case of Draco a cored or cuspy halo is expected. 

Since we shall extrapolate the dark matter mass profile well beyond the radial size of the stellar component, we 
need a description of the regime where the profile gets sensibly reshaped by tidal interactions with the dark matter 
halo of the Milky Way. We compute the tidal radius r t ;d in the impulse approximation, as appropriate for extended 
objects HHH: 



M{r m ) 



' tid 



dM 



MW 



M M w(r) dr 



M MW {r) 



,gf 



(9) 



with M(rtid) the mass of Draco within the tidal radius, and MMw( r ) the mass of the Milky Way within the galac- 
tocentric distance r; the expression on the right hand side is computed for the orbital radius of Draco r p at its latest 
pericenter passage. 

Mass models for Draco arc generated as follows: for a given functional form for the profile and for any given pairs 
of the parameters p' and a, the density profile is shifted into the form [4fj| : 



p(r) -> p(r) exp(-r/r tid ) 



(10) 



with rtid determined from Eq. @, assuming for the Milky Way a virial mass equal to 10 12 M and a NFW profile 
with concentration parameter equal to 13 J5CJ| ; r p will be taken, as a first test case, equal to 20 kpc, which is about the 
minimum pericenter radius below which tidal effects would be visible in the stellar component as well pfij , and which 
gives the most conservative estimate for the dark matter mass in Draco. We are then ready to implement Eq. and 
compare against data. 

Munoz et al.|5l| have recently made a novel compilation of l.o.s. star velocity dispersions in Draco, containing 208 
stars; they show results implementing several binning criteria, among which we resort to the one with the largest 
number of stars per bin (21 stars per bin), which is the least susce ptib le to statistical fluctuations. Using essentially 
the same data sample, but a different binning, Wilkinson et al. (52j find a sharp drop in the velocity dispersion 
corresponding to the bin at the largest circular radius, a feature that does not emerge in the analysis of Munoz et al. 
On the other hand, Lokas, Mamon & Prada question whether this sample should be further cleaned from outliers, 
i.e. stars that may not actually be gravitationally bound to Draco. We will compare separately with the data set 
from Munoz et al., i.e. in 10 bins out to a circular radius of slightly larger than 30', and the one from Wilkinson et al., 
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FIG. 3: We consider the NFW profile as reference halo model, fix the length scale parameter to 1 kpc and plot the reduced 
X 2 as a function of the density normalization parameter p . In the left panel we plot the reduced \ 2 obtained either for the 
Munoz et al. or Wilkinson et al. data sets, for a few selected values of the anisotropy parameter j3 or choosing the parameter /3 
(in the interval — 10 < j3 < 1) which at a given p' gives the smallest reduced x 2 \ other underlying assumptions (default model) 
are: radial star profile described by a Sersic profile, distance at latest pericenter passage r p — 20 kpc, mass-to-light ratio of 
3 Mq /Lq and distance of Draco d — 80 kpc. In the right panel we show that none of these latter assumptions are crucial: 
we plot as a function of p' the minimum reduced x 2 f° r P between -10 and 1, in case of radial star profile according to the 
modified Plummer model, r p increased to 80 kpc, mass-to-light ratio decreased to 1 Mq / Lq and distance of Draco varied 
within a generous range of values. 



i.e. 7 bins out to a circular radius of about 35', see Fig.0 For any mass model we consider the reduced x 2 variable: 

1 ^ (^)-<s) 2 



Xred ^ S ver y sensitive to the value of the overall normalization parameter, moderately sensitive to f3, while it is less 
sensitive to the length scale a. In Fig. [5] we show the line of sight velocity dispersion projected along the l.o.s., 
comparing to the Munoz et al. dataset and assuming the star distribution according to the Sersic profile. The best fit 
models for the three dark matter density profiles we are focusing on are set as follows: i) a NFW profile with a — 1 kpc, 
p' = 3.7 10 7 M Q kpc" 3 , r tldal = 1.7 kpc and (3 = -3.7; ii) a Burkert profile with a = 0.5 kpc, p' = 2.1 10 8 Af Q kpc" 3 , 
rudai — 2.0 kpc and (3 = —1.0; Hi) a D05 profile with a = 1 kpc, p' = 2.54 10 7 M© kpc -3 , rudai — 1-5 kpc and 
(3 = —6.3. Clearly, the dataset does not allow for a discrimination among the three models. 

In Fig. 13 we illustrate the sensitivity of the fit to some of the parameters introduced in our model, taking the 
NFW profile as reference case, and a = 1 kpc as in the best fit model: the minimum of % 2 is well defined with 
respect to p' and has a marginal shift when comparing to the data as in the binning of Wilkinson et al.; had we 
followed the suggestion of Ref. [35| to take out of the sample some of the stars that appear as outlier, the minimum 
reduced \ 2 would get below 1, but its position on the the p' axis would not change appreciably. Also shown is the 
dependence of the result upon the anisotropy parameter f3: for the NFW profile, the case of radial orbits is disfavored, 
while models with circular anisotropy give better fits. In the right panel of Fig. [31 we show instead that none of 
our additional assumptions has a significant impact on the velocity dispersion fit. In particular, there is a marginal 
effect when considering an alternative fit to the stellar profile, or when varying the assumed value for the distance of 
the Draco within the ranges of estimates quoted in the literature, or when decreasing the mass-to-light ratio of the 
stellar component to significantly smaller values. Also secondary, but slightly larger, is the effect of assuming that the 
current position of Draco is also the smallest galactocentric distance reached so far in its orbital motion, and hence it 
is the relevant radius to estimate the effects of tidal stripping (in this case, tidal radii become much larger than the 
scale radius for the stellar component). 
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FIG. 4: Left panel: for a given value of the length scale parameter a, we plot the minimum Xred obtained by varying the 
parameter p and the anisotropy parameter j3 between -10 and 1; we consider the three dark matter halo profiles introduced in 
the text, compare against the Munoz et al. dataset and refer to our default model for the other parameters. Right panel: for a 
given value of the length scale parameter a, value of p' corresponding to the minimum Xred displayed in the left panel for the 
three dark matter halo profiles. 

In Fig. 0] we show the minimum value of the reduced \ 2 > obtained taking the density normalization p' and the 
stellar anisotropy [3 as free parameters, for the three dark matter density profiles and as a function of the scale factor 
a: as clearly emerging from the Figure, the dataset does not allow for a clear discrimination in the parameter a, but 
there is, rather, a close correlation between length scale and density normalization parameter. In the right panel of 
Fig. 01 we plot the value of p' corresponding to the model with minimum \ 2 and a given scale factor a; note the huge 
span in the range of values of the logarithmic vertical scale. 

In Fig. [5] we show the tidal radii as determined assuming for the radius at the last pericenter passage 20 kpc or 
80 kpc (right panel), and values of f3 (left panel) set as in the best fit models; shallower, or less concentrated, profiles 
give equivalent fits to the data if the degree in circular anisotropy is decreased (f3 — — 10 is the minimum value we are 
scanning on; isotropic, j3 — 0, models are favored for the cored Burkert profile in the case of moderate to large values 
for the scale factor). 

B. Connections to the structure formation picture 

The possibility of discriminating among dark matter halo models increases when we take into account results from 
N-body simulations of structure formation. To make this step we need, however, to rely on a series of extrapolations. 
The first is to try to map the fit we made for a tidally disrupted object, well within the Milky Way potential well, to 
the configuration of a virialized system, unaffected by tides, of the kind described, on statistical grounds, by results of 
simulations. We refer to the prescription derived from numerical studies in |48j : let the density profile prior to tidal 
interactions be in the form: 

Pno tides(r) = p s g(r/r s ) ; (12) 
Suppose, then, that tidal interactions change it into the form: 

P(r) = /*exp(-r/r t id) p no tides (r), ( 13 ) 

assuming that the length scale parameter a in the final profile is equal to the initial scale factor r s . Comparing the 
form of Eq. (|TU1) to the one we used in the fit to the stellar velocity dispersion, i.e. Eq. (JT3J, we find p' = ftp s - The 
parameter f t is a dimensionless measure of the reduction in central density due to tidal effects; simulations indicate 
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FIG. 5: Values of /3 (left panel) and of the tidal radius (right panel, lower curves corresponding to a radius at latest pericenter 
passage equal to 20 kpc) defining the model with minimum Xred a t a given a displayed in the left panel of Fig. for the three 
dark matter density profiles. We also show in the right panel the value of the tidal radius setting the radius at latest pericenter 
passage to 80 kpc, and keeping other parameters unchanged; values for Xred change marginally, i.e. form the point of view of 
fitting the data the two cases are equivalent, but obviously all r ti d increase significantly (for reference we plot with a dotted 
horizontal line the radial size of stellar component). 



that the latter is correlated to the mass fraction of the satellite bound to the object after the effect of tides, mbnd) 
through the expression [isf : 

log ft = -0.007 + 0.35 log m b „d + 0.39 (logm b „d) 2 + 0.23 (logm b „d) 3 (14) 

(we will assume this phenomenological fitting formula to be valid for mb n d larger than about 5%). According to this 
scheme, we can uniquely assign to any best fit model with given p' and a (for an assumed pericenter radius through 
which , rtid is determined) the corresponding p s and r s , or equivalently a value for the initial virial mass of the object 
M V i r and its concentration parameter c V i r , defined as cw = R V i r /r-2- In this last step we introduced the virial 
radius R V ir, defined as the radius within which the mean density of the halo is equal to the virial overdensity A V i r 
(~ 340 at z=0) times the mean background density, and the radius r_2 where the effective logarithmic slope of the 
density profile is —2 (r_2 is equal to a for the NFW profile, 0.8a for the D05 and about 1.5a for the Burkert profile). 
In Fig. we plot M V i r for the best fit models displayed in Figs. 0] and |SJ for comparison, we also show the total 
halo mass bound to Draco after tidal stripping, and the dark matter mass within the spherical shell defined by the 
radius of the stellar component, i.e. ri um = 51'. We have referred to the two extreme choice of pericenter radii, i.e. 
20 kpc and 80 kpc; the procedure seems fairly consistent since in the two cases we get very close values for M V i r (in 
the case of a small pericenter radius and the NFW or D05 profile, at large a the fraction of mass loss becomes very 
large and extrapolations according to Eq. I|14|) becomes unreliable, so values of M V i r are not displayed). In Fig. 
we plot, for the same best fit models, the concentration parameter versus virial mass; we also show the M V i r — c V i r 
correlation as extrapolated, for the currently preferred cosmological setup j5j|, from the toy model of Bullock et al. 
54], which is tuned to reproduce the scaling found in numerical simulations for isolated halos. As far as substructures 
are concerned, concentrations are expected to be systematically larger, since substructures form, on average, in a 
denser environment with respect to isolated halos; for illustrative purposes only, we show the M v i r — c V i r scaling in 
the case of a 50% and a 100% increase in concentration. 

We have already stressed a few times that our analysis is heavily relying on extrapolations, so no firm conclusion 
can be derived; nevertheless, our results seem to indicate that we should prefer models with an intermediate M vir , 
say 10 9 Mq , corresponding to a of the order of 1 kpc for the NFW and D05 profiles and about 0.5 kpc for the Burkert 
profile (such cases are those that we have been chosen as reference models in Fig. 0) and that the range of length scale 
values allowed in Figs. 0] and [S] is probably a very generous one, with values at the lower and upper ends which should 
be most likely dropped. The range of models we are indicating here as preferred by the NFW profile is analogous 
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FIG. 6: We display, for the minimum Xred models illustrated in Figs. 2] and values of the dark matter halo mass within the 
radial size of stellar component (lower curves) and of the total mass in the dark matter component for the calculated tidal 
radius (medium curves), assuming a pericenter radius of 20 kpc (left panel) and 80 kpc (right panel). We have also performed 
an extrapolation to estimate the initial virial mass of Draco, i.e. the mass associated to it before sinking deep into the potential 
well of the Milky Way and the loss of a large fraction of such initial mass due to tidal effects. 
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FIG. 7: We display, for the minimum Xred models illustrated in Figs. 2] and the extrapolated values of the virial mass and 
concentration parameter according to the prescription for the response to tidal interactions introduced in |4gj| and fitted to 
numerical simulations. Also shown is the extrapolation with the Bullock et al. prescription of the correlation M v ir — Cvir for 
isolated halos, and assuming a 50% or a 100% increase in concentration in the case of subhalos. 



to the one suggested in Ref. [25j, although the two approaches differ. In particular we will not implement here a 
constraint from the age of Draco stellar population which is used as guideline in |25j| : to do that we would need to 
build a subhalo mass function for the Milky Way matching the observed satellite pattern, and to model star formation 
within subhalos, two steps which are not very well understood and on which the degree of extrapolation would be 
inevitably much more drastic than what we have accepted so far. 
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III. THE GAMMA-RAY SIGNAL FROM WIMP ANNIHILATIONS IN DRACO 



WIMPs have a small but finite probability to annihilate in pairs, giving rise to potentially observable standard 
model yields. Two ingredients intervene in fixing source functions: on the one hand, the annihilation cross section, 
branching ratios and spectral distributions for the yields are specified in any given particle physics scenario embedding 
the WIMP; on the other hand, source functions scale with the number density of WIMP pairs, i.e. in the case we 
are considering here, they are proportional to the square of the dark matter mass density in Draco. Since photons in 
the energy range we are interested to, i.e smaller than few TeV, propagate through the interstellar medium without 
being absorbed, predictions for the induced gamma-ray fluxes are straightforward and simply involve an integral of 
the source along the line of sight; the expression for the flux per unit energy and solid angle, is usually casted in the 
form: 



where (cry) is the WIMP annihilation rate at zero temperature, M x the WIMP mass and the sum is over all kine- 
matically allowed annihilation final states /, each with a branching ratio Bf. It is beyond the scope of the present 
analysis to review the ranges of values and the model-dependent determination of these parameters, as well as of the 
gamma-ray spectral distributions dN^/dE, topics which have been vastly discussed in the recent literature; we will 
mainly refer here either to a toy- model in which we pick particular values for M x and (cry), and assume to have only 
one dominant annihilation channel: it is useful to consider the case of a soft annihilation channel such as b — b pair, 
and to contrast it with the hard r + — t~ final state. As we showed in these toy models are well justified in the 
context of solidly motivated theoretical grounds, for instance within the paradigm of neutralino dark matter. For 
definiteness, and for illustrative purposes, we shall also make use of special, well-studied, benchmark supersymmetric 
models, as we did for the case of our analysis of the multi- wavelength emissions from Coma in Ref. Q] . 
In Eq. (|15Jl the dependence on the halo profile has been factorized out defining: 



where O is the direction of observation and the average is over the angular acceptance (or the angular resolution) 
of the detector AfL In Fig. |H1 we plot the range of the expected values for J towards the center of Draco, for two 
sample values of Af2 and within the minimum x 2 halo models selected in the previous Section. Confirming other 
recent analyses |l3l l55l 156$ . our results show that there is a very small spread in the prediction for J, even referring 
to significantly different dark matter halo shapes and even for small angular acceptances: within a factor of few and 
in units of GeV 2 cm" 6 kpc, J is about 100 for A£l = 10~ 5 sr and about 1 for A£l = 10~ 3 sr. Such small spread is 
in contrast to what one finds in the analogous estimate when considering the Milky Way galactic center as a source 
of gamma rays from dark matter annihilation. One can apply the same procedure of fitting different halo profiles to 
the Milky Way dynamical constraints and then extrapolate their radial scalings down to the innermost parsec or so; 
the focus is on the eventual sharp dark matter density enhancement which could be present in the Galactic center 
region: for singular profiles the values of J one derives may be very large, up to about 10 4 -10 5 for NFW profiles 
and AQ = 10~ 5 sr (see, e.g. ^j|; note however that in Ref. 01 a dimensionless J is adopted and to translate values 
quoted therein into those for the definition adopted here, one should scale them by the factor 0.765 GeV 2 cm~ 6 kpc), 
but drop dramatically, with a decrease as large as four orders of magnitude, when considering less singular or cored 
profiles, in the case of Draco, the distance of the source is much larger and the l.o.s. integral involves an average over 
a large volume, smoothing out the effect of a singularity in the density profile; at the same time, however, the mean 
dark matter density is on average fairly large for any profile, since the dark halo concentration parameter is large. 

In its all-sky survey, EGRET has accumulated a limited exposure towards Draco. A report on the collected data 
is given in [l||; the analysis aims at the identification of a point source at the center of Draco; seven energy bins are 
considered, each with the appropriate angular cuts, no point source is found, and the photon counts are consistent with 
the expected flux from diffuse emission, except for a 2 event "excess" in the energy range between 1 and 10 GeV, with 
a total of 6 events found versus 4.1 expected in standard background scenario (notice that no statistical evidence for 
such "excess" is claimed in |l^ or in the present analysis). In Fig.[!|]we show, for a given WIMP mass, the value of the 
annihilation cross section required for a flux matching the 2 events in EGRET between 1 and 10 GeV, for exposures 
and angular cuts as specified in the data analysis, assuming our reference NFW best-fit halo model and a b — b (left 
panel) or r + — r _ (right panel) annihilation channels. Also shown in the Figure are expected sensitivity curves with 
GLAST, the next gamma-ray telescope in space, and for upcoming observations of Draco with ground-based ACT 
telescopes. 




(15) 




(16) 
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FIG. 8: We display, for the minimum Xred models illustrated in Figs.^JandQj] the line of sight integral function J, towards the 
center of Draco and for two different angular acceptance AS1. 




FIG. 9: We show the WIMP mass - WIMP pair annihilation cross section plane, assuming a b — b (left) and r + — r~ 
(right) dominant annihilation final state. The solid black line indicates the value of the annihilation cross section required 
for a flux matching the 2 events in EGRET between 1 and 10 GeV, for exposures and angular cuts as specified in the data 
analysis, assuming our reference NFW best-fit halo model; The black dotted and dot-dashed lines show the limits from the 
flux of, respectively, antiprotons and positrons in the Milky Way halo (the two lines corresponding to two different cosmic ray 
propagation setups); The solid red, dark blue and light blue curves indicate the projected sensitivities of GLAST, MAGIC and 
VERITAS (see the text for more details). 



Regarding the GLAST detector, we refer to an updated simulation of the instrument performance [57j: we refer to 
the energy dependent sensitivities of the two LAT sections, the thin (or front) section of the tracker (peak effective 
area above 1 GeV of about 5500 cm 2 , 68% containment angle varying between 0.6 deg at 1 GeV and 0.04 deg at 
100 GeV) and of the thick (or back) section of the tracker (peak effective area above 1 GeV of about 4500 cm 2 , 68% 
containment angle varying between 1 deg at 1 GeV and 0.07 deg at 100 GeV). To estimate the background, we include 
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an extragalactic component at the level found in EGRET data [5g, extrapolated to higher energy with a E~ 2A power 
law, plus a galactic component scaling like E~ 2 - 7 (such scaling is expected from the decay of it generated by the 
interaction of primary protons with the interstellar medium; we are neglecting an eventual IC component, since, if 
present, such term is most likely already included as misidentified extragalactic) and normalized in such way that, 
together with the extragalactic component, it gives the 6 events above 1 GeV detected by EGRET (assuming just for 
4.1 events for the background level does not change significantly our projected sensitivities). We consider a 5 years 
exposure time, in an all-sky survey mode for which the effective area in the direction of Draco is, on average, about 
30% of the peak effective area (area when the source is at the zenith of the instrument). Finally, we define a x 2 
variable as: 



nbi» a N 2 

y 2 = V J (17) 

x ^ Ns; + Nbj K ' 

3=1 J J 

where Nsj and Nbj stand for the number of signal and background events in each bin, restricting to bins with more 
than 5 signal events. The bin selection should in principle be optimized model by model; in general we find that 
it is a good choice to take three bins per decade in energy (two or more bins are grouped into one in case this 
procedure gives 5 signal events; this sometimes happens in the highest energy bin included in the sum above), while 
at any given energy we integrate over a solid angle which is the largest between the PSF set by the 68% containment 
angle (full energy dependence included for each section of the tracker) and the solid angle which maximizes the ratio 
^ s /y/Wl + $h. Sensitivity curves are given in Fig. as 3 a discovery limits; the latter are found to be, with the 
present accurate modeling of the detector, slightly less promising than the analogous curves obtained in other recent 
estimates [55ll56| . 

Regarding Air Cherenkov Telescopes (ACTs), we consider the detection prospects with instruments in the northern 
hemisphere, i.e. MAGIC which is currently taking data, and VERITAS, which will be completed soon extending 
the current single mirror telescope to an array of at first four, then later seven telescopes [(hJ . We assess the discovery 
sensitivity of the two ACTs using a low energy threshold of 100 GeV, and the effective collection area as a function of 
energy recently quoted by the two collaboration in Ref. 59, 60] ; the main sources of background for ACTs correspond 
to misidentified gamma-like hadronic showers and cosmic-ray electrons. The diffuse gamma-ray background gives a 
subdominant contribution to the background, which we also took into account using the same figures outlined above 
for the space-based telescopes background. We use the following estimates for the ACT cosmic ray background 61] : 



-2.3 



~^-(E > E ) = 6.1 x 10- 3 e had Iy^v) cnrV 1 !*- 1 , (18) 



^(E>£b)-3.0X10-( T ^)' cm-W, (19) 

where Chad parameterizes the efficiency of hadronic rejection, which we assume to be at the level of 10% [sol Ifffll]. As 
above, we proceed with an optimized binning of the energy range of interest (extending from the energy threshold up 
to the WIMP mass), compute the number of signal and background events in each bin, and require that the resulting 
X 2 (evaluated according to the analogue of Eq. I|17|l) gives a statistical excess over background. 

The models for which we predict a detectable flux have fairly large cross sections, still compatible but in the high 
end of models with a thermal relic density, as computed in a standard cosmological scenario, which matches the 
observed dark matter density in the Universe [H^|. see, e.g., Q (another possibility is that we refer to models with 
non-thermal relic components, such as from the decay of moduli fields, or to modified cosmological setups affecting the 
Hubble parameter at the time of WIMP decoupling, see, e.g., Large annihilation cross sections give enhanced 

signals for any indirect dark detection technique, in particular we need to check whether this picture is compatible 
with the antimatter fluxes measured at Earth: in fact, pair annihilation of WIMP in the halo of the Milky Way is 
acting as a source of primary positrons and antiprotons which diffuse in the magnetic field of the Galaxy, building 
up into exotic antimatter populations. Current measurements of the local antiproton flux are consistent with the 
standard picture of antiprotons being secondary particles generated by the primary cosmic ray protons in spallation 
processes [6^ ]: on the other hand, a weak evidence of an excess in the positron flux has been claimed |64l l65l] . in a 
picture that is going to become increasingly clearer with the on-going measurements in space by the recently launched 
Pamela detector [661 ] . We estimate the induced flux of positrons and antiprotons (no antiprotons are generated in 
the t + — t~ channel), referring to the same Milky Way halo model we have introduced to estimate tidal effects on 
Draco, and to the diffusive convective model for the propagation of charged particles implemented in the DarkSUSY 
package [(JlJ ■ Parameters in the propagation model are chosen in analogy with a standard setup |6^ in the GALPROP 
propagation package or the most conservative choice suggested in Ref. [7(j which can still reproduce ratios of 
secondaries to primaries as measured in cosmic ray data while minimizing the flux induced by WIMP annihilations: 
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FIG. 10: As in Fig. 03 but keeping the WIMP mass fixes at 100 GeV, and varying the halo scale length parameter a for the 
three dark matter density profiles considered here. 



these give, respectively, the lower and upper curves plotted in Fig. and corresponding to the 3 a limits on the 
annihilation cross section obtained by comparing the WIMP-induced fluxes to a full compilation of present data on 
the local antiproton and positron fluxes. The values displayed should should not be taken as strict exclusion curves, 
since we are not doing a full modeling of the uncertainties in the propagation model, nor scanning on more general 
configurations of the Milky Way dark matter halo; relaxing them by a factor of 2 to 5 or maybe even larger should be 
relatively straightforward; they can be however taken as a guideline to see that models in such region of the parameter 
space should be testable with higher precision antimatter data, while models at the EGRET level are most probably 
already excluded by current data. 

Finally in Fig. E3 w e show how flux levels and projected sensitivities scale scanning over our sample of halo models 
for Draco, for a sample WIMP mass and the b — b final state; antiproton limit levels are also plotted for comparison. 

This work which is based on realistic estimates of both the dark matter density distribu- 
tion in Draco and of th e WIMP physical set up leaves the early Cactus claims (see, e.g., 
http://ucdcms.ucdavis.edu/solar2/results/Chertok.PANIC05.pdf) of a detection of gamma-ray emission from 
Draco apart. 



IV. MULTI- WAVELENGTH SIGNALS FROM DRACO 



The following step is to extend our analysis to the radiation emitted at lower frequencies. For this purpose, we 
need to track the injection of electrons and positrons from WIMP annihilations in Draco, as well as their propagation 
in space and energy; it will then be possible to make predictions for the induced synchrotron and Inverse Compton 
radiations. Our starting point is the assumption that, in analogy to more massive objects such as the Milky Way itself, 
there is a random component of interstellar magnetic fields associated with Draco and that it is a fair approximation 
to model the propagation of charged particles as a diffusive process. In this limit we can calculate the electron and 
positron number densities implementing to the following transport equation: 



d dn e 
di~d~E 



= V 



D{E,x)V 



dn e 
~dE 



d_ 

dE 



b(E,x 



dn e 
"dE 



Q e {E,x) , 



(20) 



where Q e is the electron or positron source function from WIMP annihilations: 

dNf 



Qe(E,x) 



—T 



2 A/2 ^ dE 

x f 



(E)B f p\x), 



(21) 



14 




FIG. 11: Left panel: mean diffusion distance V Av covered by an electron while losing energy from its energy at emission E' 
to its energy at interaction E; three different E' are considered, as well as a conservative choice for propagation parameters 
(labeled by "set # 1" in the plot) or a more extreme choice (labeled by "set # 2" in the plot), see the text for details. Right 
panel: we plot, as function of the distance from the center of Draco r, the integral over volume up the radius r of the electron 
number density dn e /dE(r, E), for three values of the energy E and for the two choices of propagation parameters already 
considered in the left panel; also shown is the case when dn e /dE is computed assuming negligible spatial diffusion. For each 
energy, integrals over volume are normalized to the integral over volume up to the assumed radius of diffusion for Draco of 
dn e /dE in the limit in which spatial diffusion is neglected; we have chosen the reference NFW model for the dark halo, and 
a WIMP of mass 100 GeV annihilating into bb. As it can be seen, compared to the case when spatial diffusion is neglected, 
in the actual cases applying to Draco there is a sharp deficit of electrons in the inner region of even within the total diffusion 
volume. 



while D is the diffusion coefficient and b the energy loss term. In Ref. [l| we derived the analytic solution to this 
equation in case of a spherical symmetric system and for D and b that do not depend on the spatial coordinates. 
We refer here to a time-independent source and consider the limit for an electron number density that has already 
reached equilibrium; the solution takes the form: 

dn 1 f Mx - 

-di (r ' E)= b(E)JE dE ' G ^ v - v ')Q^ E ') ( 22 ) 

with: 

[47r(A«)]V2 n £^J J Q r n [ l \ AAv J ^\ AAv J J p 2 (r) 

In the Green function G the energy dependence has been hidden in two subsequent changes of variable v = 
duD{u) and u = /J' max if(E') > ^ ne ra( iial integral extends up to the radius of the diffusion zone r/j at which 
a free escape boundary condition is imposed, as enforced by the sum over n, having defined r„ = (— l) n r + 2nrh- 
Whenever the scale of mean diffusion V Av, covered by an electron while losing energy from energy at the source E' 
to the energy when interacting E, is much smaller than the scale over which p 2 has a significant variation, G is close 
to 1 and spatial diffusion can be neglected; in Ref. it was shown that this limit applies in the case of the Coma 
cluster. For dSph we find that, most likely, we are in the opposite regime. 

To model electron and positron energy losses we choose as a reference value for the magnetic field = 1 /iG and 
a thermal electron density of 10~ 6 cm -3 : these values are derived from radio observations of dwarf galaxies similar 
to Draco at 5 GHz (Klein et al. 1992) and from the assumption that the ROSAT PSPC X-ray upper limit on Draco 
(Zang & Meurs 2001) is due to thermal bremsstrahlung, respectively. For the diffusion coefficient we assume the 
Kolmogorov form D(E) = D Q /Bl /3 (E/l GeV) 1/3 , fixing the constant D Q = 3- 10 28 cm 2 s 1 in analogy with its value 
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FIG. 12: Radio flux density spectrum (left panel) and surface brightness distribution at the frequency v — 4.9 GHz for a sample 
WIMP model of mass 100 GeV annihilating into bb with a cross section tuned at the level to give 2 events in the EGRET 
gamma-ray telescope. Results are given for the two choices of propagation parameters already considered in Ref. 1111 and for 
the surface brightness, for illustrative purposes only, in case when spatial diffusion is neglected. Surface brightness are plotted 
for a 3 arcsec angular acceptance, corresponding to the VLA angular resolution at the time a search with this instrument for 
a point source at the center of Draco was performed (the obtained upper limit is plotted in the figure) . 



for the Milky Way; finally our guess for the dimension of the diffusion zone is that, again consistently with the picture 
relative to the Milky Way, is about twice the radial size of the luminous component, i.e., here, 102 arcmin (we will 
refer to this set of propagation parameters as set #1). In the left panel of Fig.^Jwe show that, in this setup, electrons 
and positrons lose a moderate fraction of their energies on scales \J Av that are comparable to the size of the diffusion 
region, i.e. spatial diffusion as a large effect. Even referring to an extreme model in which the diffusion coefficient 
is decreased of two orders of magnitudes down to D = 3 • 10 26 cm 2 s _1 (this would imply a much smaller scale of 
uniformity for the magnetic field), adding on top of that a steeper scaling in energy, D(E) = Do (E/l GeV) 0,6 (this 
is the form sometime assumed for the Milky Way; we label this propagation parameter configuration set #2), scales 
V At; are decreased but remain still relatively large. In the right panel of Fig. ^] we consider a WIMP model with 
mass 100 GeV annihilating in the bb final state within our reference NFW halo model for Draco. We show integrals 
over volume within the radial coordinate r of the equilibrium number density dn e /dE, for a few values of the energy 
E, and for the set of propagation parameters #1 and #2, as well as the results corresponding to the assumption 
that spatial diffusion can be neglected. All integrals are normalized to the integrals over the whole diffusion region 
of dn e /dE for the corresponding energy E and assuming negligible spatial diffusion: we deduce from the figure that 
for set #1 there is a depletion of the electron/positron populations with a significant fraction leaving the diffusion 
region, while for set #2 they are more efficiently confined within the diffusion region but still significantly misplaced 
with respect to the emission region. 

For a given electron/positron equilibrium distribution we can infer the induced synchrotron and inverse Compton 
emissions. In the limit of frequency v of the emitted photons much larger than the non-relativistic gyro-frequency 
= eB / {2-Kmc) ~ 2.8-B^ Hz, the spontaneously emitted synchrotron power takes in the form [7l|: 

^synch {v, E, r) = [ d9^^-2irVir mciy sin 9 F {x/ sm9) , (24) 
Jo 2 

where we have introduced the classical electron radius tq — e 2 /(mc 2 ) — 2.82 • 10 -13 cm, and we have defined the 
quantities x and F as: 



2v 

x : 



3i"o7 2 



7*V 2 



3/2 

(25) 
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and 



/DO 
dzK 5/s (z) ~ 1.25i 1/3 exp (-i) [648 + t 2 ] 1/12 . (26) 

Folding the synchrotron power with the spectral distribution of the equilibrium number density of electrons and 
positrons, we find the local emissivity at the frequency v. 



, Mx ( dn„- dn„+ , 
Jsynch {y, r) = / dE ( + ) P synch (u, E, r) . (27) 



In Fig.^Jwe consider a reference WIMP model with a mass of 100 GeV, annihilating into b—b with a cross section 
at the level to induce a gamma-ray flux matching the 2 events in EGRET between 1 and 10 GeV for a dark matter 
distribution as in our reference NFW model. In the left panel we plot the radio flux density spectrum integrated over 
the whole diffusion volume: 

S synch (v) = /^ jsy ^' r) , (28) 

with d the distance of Draco. The spectrum is significantly flatter for the propagation parameter set #1, since the 
peak in emitted photon frequency is linearly proportional to the energy of the emitting particle, and electrons and 
positrons tend to escape from the diffusion box while losing energy, rather than remaining confined within it and 
giving a signal which piles up at lower frequencies. We then introduce the azimuthally averaged surface brightness 
distribution: 

I synch { V ,e,AQ)= [ <Kl [ dl Jsync ^ Uj) , (29) 

JAU Jl.o.s. 4 7T 

where the integral is performed along the l.o.s. /, within a cone of angular acceptance AQ centered in a direction 
forming an angle O with the direction of the center of Draco. In the right panel of Fig. 1121 we plot surface brightness 
integrated over a cone of 3 arcsec width, corresponding to the tiny angular acceptance of the VLA at the time it was 
used to perform a searches for point radio sources in the central 4 arcmin of Draco [2£j ; no source was found and the 
corresponding upper limit is plotted in the figure. To illustrate how the shape of the signal changes compared to that 
of the source function, we also plot the surface brightness which we would obtain in the limit of no spatial diffusion. 
In Ref. radio fluxes are computed in this limit and the VLA measurement is exploited to exclude WIMP models; 
we have demonstrated in our discussion that the limit of no spatial diffusion is not likely to hold in case of Draco and 
the figure illustrates the fact that, with present data, limits on the model stemming from radio frequencies are less 
constraining than in the gamma-ray band. The figure illustrates also another point: the gamma-ray flux retraces the 
WIMP annihilation source function; in the example we have considered, even with the angular resolution at which 
future observations will be carried out, Draco would appear as a single point source. On the other hand, in the radio 
band the signal is spread out over a large angular size, standing clearly as diffuse emission. 

No search for a diffuse radio emission from Draco has been performed so far. Even with a future next generation 
radio telescope the quoted sensitivities do not apply for an extended source. To address the discovery potential of 
such apparata for the signal we are analyzing, we make a simple extrapolation on the quoted point source sensitivity 
I p {v) for a reference angular resolution AQ p and integration time At p , assuming a homogeneous background: 



I mm (y,£n)=I p (v) -A—= ^j=l . (30) 
w Ip v At 

Reference values for EVLA in phase I [73|. at v — 5 GHZ are I p {v) = 0.8 /LtJy for At p = 12 hr and AQ p — 0.4 arcsec; 
for LOFAR observations are al lower frequency up to v = 200 MHZ, for which I p (v) = 0.03 mJy with At p — 1 hr and 
AJ7 p = 0.64 arcsec |74j. To infer the projected sensitivity limit, for each WIMP model, halo profile and frequency of 
observation, we compute the value for the angular acceptance AQ at which / S ynch(^, AS1) / \J AS1 is maximal; we also 
assume as exposure time At = 8 hr. In Fig. 1131 we show results for the sensitivity curves in the plane annihilation cross 
section versus WIMP mass, for our reference NFW profile, the conservative set #1 for propagation parameters (set 
#2 would give much more favorable results) and a value of the magnetic field equal to 1 /iG. The figure shows that, for 
this choice of parameters, we predict that WIMP models that are at the level of being detected by GLAST or ACTs in 
gamma-rays, should also give a detectable radio flux, possibly with an even higher sensitivity in favorable propagation 
configuration. However, this last conclusion is more model dependent, since some of the parameters are crucial for 
the estimate of the radio flux. In Fig. ^] we illustrate this point, fixing the WIMP mass to 100 GeV, varying instead 
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FIG. 13: The projected sensitivity of future diffuse radio source searches from the direction of Draco, in the same setup as 
m Fig. M For comparison, we also indicate the sensitivity of gamma-ray search experiments, and the constraints from the 
antiproton and positron fluxes (the two lines corresponding to two different propagation setups for the Milky Way). 
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FIG. 14: The projected sensitivities of future searches for a diffuse radio source as a function of the magnetic field in Draco, 
in the same setting as in Fig. |5] with a WIMP mass set to 100 GeV. For comparison, we also show the projected sensitivity of 
GLAST. 



the value of magnetic field in Draco: the trend is obviously that the larger the magnetic field, the higher the induced 
radio flux, but the dependence is not trivial since the magnetic field enters both in the propagation of electrons and 
positrons, and in the emission of synchrotron radiation. Finally in Fig. 1151 we examine the dependence of the radio 
sensitivity curves on the model describing the dark matter halo in Draco: we find scalings that are analogous, although 
different in fine details, to those sketched previously for the gamma-ray fluxes and the corresponding sensitivity of 
the GLAST satellite. 

The inverse Compton emission on the cosmic microwave background and on starlight fills the gap in frequency 
between radio and gamma-ray frequencies. Let E = jrrieC 2 be the energy of electrons and positrons, e the target 
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FIG. 15: LOFAR projected sensitivities, for a given WIMP model and value of the magnetic field in Draco, as a function of 
the length scale parameter a in the class of halo models selected in Sec. [HI compared to the GLAST projected sensitivity. 



photon energy and E 1 the energy of the scattered photon; the Inverse Compton power is obtained by folding the 
differential number density of target photons with the IC scattering cross section: 

P lc {E 1 ,E) = cE 1 J den{e)(j{E 1 ,e,E) (31) 

where a{E 1 ,e,E) is the Klein-Nishina cross section |7l| and n(e) is the differential energy spectrum of the target 
photons; for simplicity we will assume that the starlight spectrum has the shape of a black body with temperature 
T = 0.3 eV Such value of the temperature has been estimated on the basis of the fact that the major part of the 
Draco star are halo stars somewhat below the turnoff point of the subdwarf main sequence (Odenkirchen et al. 2001). 
The effective temperature in the HR diagram with Fc/H=-2.0 dex and with an age of ~ 12 Gyr is of the order of 
T rj 3300 — 3500 K, which is equivalent to an energy of « 0.28 — 0.3 eV. Folding the IC power with the spectral 
distribution of the equilibrium number density of electrons and positrons, we get the local emissivity of IC photons 
of energy E~ ( : 

he (E~ n r) = JdE + ^) Pic {E n ,E) (32) 

and the azimuthally averaged surface brightness distribution: 

7i C (*/,e,Afi) = / dtt f dl jlC y " . (33) 

In Fig. ^3 we plot the sample multi- frequency seed of the emission in Draco due to WIMP annihilations, imple- 
menting our reference NFW halo model and reference values for the WIMP mass, the magnetic field and the various 
propagation parameters. The WIMP pair annihilation rate has been tuned to give a gamma-ray signal at the level 
of the EGRET measured flux; the displayed surface brightness is in the direction of the center of Draco and for an 
angular acceptance equal to the EGRET angular resolution, i.e. not optimized for future observations (we should 
have in fact considered different solid angles at different wavelengths). As apparent, there is a significant component 
in the X-ray band due to inverse Compton on the microwave background radiation, while the contribution on starlight 
is essentially negligible. Scaling of signals with the assumed value of the magnetic field are also displayed in the right 
panel. 

While so far we have considered simple toy models for the WIMP accounting for the dark matter halo in Draco, 
in Fig. 1171 we consider a few explicit realizations of this scenario within the constrained minimal supersymmetric 
extension to the Standard Model (cMSSM), picking among the models studied in [7{| those better exemplifying the 
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FIG. 16: Detailed multi- wavelength spectrum for a 100 GeV WIMP annihilating into b — b (left), and the effect of varying 
the magnetic field strength. The WIMP pair annihilation rate has been tuned to give a gamma-ray signal at the level of the 
EGRET measured flux. 



widest range of possibilities within that particular theoretical setup. All the models are fully consistent with accelerator 
and other phenomenological constraints, and give a neutralino thermal relic abundance exactly matching the central 
cosmologically observed value j^. We adjusted here the values of the universal soft supersymmetry breaking scalar 
mass mo given in |7^| in order to fulfill this latter requirement, making use of the latest Isajet v. 7. 72 release and of 
the DarkSUSY package [S^. The values of the cMSSM input parameters for the various models are given in Tab.|U 
(see also Ref. 0]). Each benchmark model correspond to a different mechanism responsible for the suppression of the 
otherwise too large bino relic abundance: B' lies in the bulk region of small supersymmetry breaking masses, and 
gives a dominant b — b final state; D' corresponds to the coannihilation region, and features a large branching ratio 
for neutralino pair annihilations in r + — t~ ; E' belongs to the focus point region, with a dominant W + — W~ final 
state, and, finally, K' is set to be in the funnel region where neutralinos rapidly annihilate through s-channel heavy 
Higgses exchanges, dominantly producing b — b pairs as outcome of annihilations. Not unlike what we found in the 
case of the multi- wavelength analysis of neutralino annihilations in the Coma cluster (see fig. 25 in Ref. Q), the most 
promising among the four benchmark models of Tab.[I]is model E', featuring a large pair annihilation cross section to 
begin with; the less promising model is instead model D', for which the mechanism suppressing the neutralino relic 
abundance in the Early Universe, stau coannihilations, is not associated to pair annihilations of neutralinos today. 

Lastly, we found that the SZ effect produced by DM annihilation in Draco, even though is a definite probe of 
the DM annihilation in such cosmic structures (see, e.g., |9j,[7j|) is quite low when we take into account the spatial 
diffusion of secondary electrons: we find, in fact, that the SZ signal towards the center of Draco is negligible even 
when we normalize the gamma-ray signal at the level of the EGRET upper limit. 
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TABLE I: The input parameters of the four cMSSM benchmark models we consider here. The units for the mass parameters 
are GeV, and the universal trilinear coupling A a is set to for all models (see 0,|7^] for details). 
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FIG. 17: Multi- wavelength spectra for the four benchmark models described in the text. The prediction is shown for the best 
fit NFW profile, and for a mean magnetic field equal to 1/xG. 

V. A BLACK HOLE AT THE CENTER OF DRACO? 

There is one further effect which could change substantially our picture: if a black hole is present at the center 
of Draco, and had it formed through an adiabatic accretion process, the ambient dark matter population would 
have experienced a sharp increase in its density profile, turning into a "spike" of dark matter, with a dramatic 
enhancement in the dark matter annihilation rate at the center of Draco. Such spike was originally proposed for 
the Milky Way [H EE El E3 

in connection its central black hole, which has a mass of about 3 10 6 M Q , and, more 
recently, it has been extrapolated to small mass halos 0, |i3 > including substructures within the Milky Way dark 
matter halo, eventually embedding black holes of intermediate mass, in the range between 10 2 to 10 6 M Q . 

There is a strong observational evidence for the existence of super-massive black holes (in mass range between 10 6 
and 1O 9 M ), without however a detailed understanding on how those objects form, or on the mechanism enforcing 
the observed correlations with properties of the hosting halos. In one of the proposed scenarios, these two issues 
are addressed in terms of pre-existing intermediate-mass black hole seeds, forming in turn in proto-galaxy environ- 
ments [13, El : a significant population of these smaller mass objects would still be present in Galaxy-size halos, 
most likely associated to substructures which have not been tidally disrupted, while merging into the halo. Their 
presence in the Milky Way halo would be very hard to prove in terms standard astrophysical observations. In par- 
ticular, there is no evidence for the presence of a black hole at the center of Draco: it is reasonable to expect that 
a black hole, being in such a gas poor environment, would be in a dormant phase, rather than in an accreting and 
luminous one. In Fig. ^| we sketch the dynamical response of adding a black hole of given mass on top of the mass 
models introduced in Section HJ (the response of the dark matter profile, as specified below, is included): the fit of 
the star velocity dispersion is not sensitive to black holes of mass smaller than about 10 5 M Q , slightly improves for 
masses around a few times 10 6 M Q , while the presence of black hole of mass larger than about 10 7 M Q is dynamically 
excluded. 

We will take a phenomenological approach and make the hypothesis that a black hole of given mass Mbh has 
formed adiabatically at the center of Draco. The process turns an initial (i.e. before the black hole has accreted the 
bulk of its mass) dark-matter density profile scaling as /Oj(r,) oc r i 7 into a final profile of the form pf(r/) oc rj : in 
a simplified system with all dark matter particles on circular orbits, it is easy to show that conservation of mass and 
angular momentum imply that A = (9 — 2j) /(4 — 7) jz3, EE EH 1 i- e - that the final profile is significantly steeper than 
the initial; this results holds also in a general setup. To derive the right normalization, on the other hand, one has to 
refer to the full phase space distribution function for the dark matter profile and implement the appropriate adiabatic 
invariants. We refer here to the procedure outlined in |7S|: in the same paper it is shown that, since the growth of the 
spike depends on the existence of a very large population of cold particles at the center of the dark matter system, 
where the black hole is adiabatically growing, large spikes form for singular profiles, which embed such large number 
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FIG. 18: Reduced \ 2 f° r a fit of the star velocity dispersion data of Draco in the binning of Munoz et al. under the hypothesis 
that a black hole of given mass is present at the center of Draco; the dark matter profiles are described by a NFW or a D05 
profile with scale factor a = 1 kpc and, respectively, p' equal to 10 7 , 2 ■ 10 r , 3 • 10 7 , 3.72 • 1O 7 M0 kpc -3 (from top to bottom in 
the figure; the last value corresponds to best fit in the case without the black hole) and 5 • 10 6 , 10 7 , 2 ■ 10 7 , 2.54 ■ 10 7 M© kpc -3 
(again, from top to bottom in the figure, with last value being the best fit in the case without the black hole). 



of cold statesr, while it does not for cored profiles for which it is not the case. We will discuss then only the case for 
the NFW profile and the D05 profile. 

In Fig. ll9l we plot the line-of-sight integral function J, we have introduced in Eq.^as relevant quantity for predictions 
on the gamma-ray flux, as a function of the black hole mass and for a given value of the WIMP annihilation cross section 
av, or vice versa. The value of av enters critically since the very singular spike density profile has to be extrapolated 
down to the radius at which a maximal WIMP density is enforced by WIMP pair annihilations, i.e. refGondolo:1999ef: 



M x 
(av)(to -t f ) 

where to is the present time and tf the formation time of Draco. As it can be see, enhancements in the gamma-ray 
flux of even four orders of magnitude are at hand; scalings in black hole mass and av are analogous for the two halo 
models considered here. 

The spike is confined in a very tiny portion of the halo, essentially the region within which the black hole dominates 
the potential well of the final configuration, i.e. smaller than 1 pc even for the heaviest black holes we are considering. 
The induced gamma-ray source would appear as a point source even at future telescopes with improved angular 
resolution. On the other hand, analogously to the effect we have already discussed for the standard dark matter halo 
component, the emitted electrons and positrons diffuse out of the central region and give rise to radio and Inverse 
Compton signals on a very wide angular size. In Fig.[20|we consider, for a few sample masses for the central black hole 
and one reference annihilation cross section, the induced radio surface brightness in the same configuration displayed 
in Fig. 1121 (propagation parameters in set #1). 

In Fig. 1211 we show the scaling of future expected sensitivities in the plane annihilation rate versus black hole mass, 
for two sample WIMP masses, the 6—6 annihilation channel, the reference NFW halo profile and set of propagation 
parameters. As already mentioned above, the effect of the adiabatic black hole growth is more dramatic for WIMP 
models with smaller annihilation cross section; for comparable annihilation cross section the enhancement in the 
signal is larger at radio wavelengths than for gamma-rays. There is a 1/4-7T mismatch since we are essentially adding 
a point source at the center of the system: we detect only the gamma-rays emitted in our direction, while all emitted 
electrons and positrons pile up into the population giving rise to the radiation at lower frequencies. 
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FIG. 19: Left panel: integral of the square of dark matter density along the line of sight towards the center of Draco and 
averaged over an angular acceptance of 10~ 5 sr, in case of a central enhancement of the dark matter density due the adiabatic 
formation at the center of the system of a black hole of given mass; the background dark matter profile are those considered 
in the left panel and for p' equal to 3.72 • 1O 7 M kpc" 3 (NFW profile) and 2.54 • 10 7 M o kpc" 3 (D05 profile). The value of 
the WIMP pair annihilation cross section and mass enter in the estimate of central maximal WIMP number density set by 
pair annihilations. Right panel: same as in the right panel but for a few values of the black hole mass and varying the pair 
annihilation cross section. 




FIG. 20: We show the radio (at v — 4.9 GHz) surface brightness of Draco in the case of an adiabatic growth of a black hole 
of given mass Mbh = 10 2 , 10 4 and 10 6 Mq, as labelled. A reference neutralino bb model with M x — 100 GeV and av = 10 -25 
cm 3 s _1 with a magnetic field _B M = 1 fiG are adopted here. 



VI. CONCLUSIONS 



Following a program of thorough investigation of the multi-wavelength yields of WIMP dark matter annihilations 
started with the case of the Coma cluster in Ref. [lj, m this paper we analyzed the case of the nearby dSph Draco. 
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FIG. 21: We show the projected sensitivity of radio (LOFAR, EVLA) and gamma-ray (EGRET, GLAST) experiments observing 
Draco as a function of the black hole mass Mbh for two choices of the neutralino mass: M x = 100 GeV (left panel) and M x = 500 
GeV (right panel). A reference case of a NFW dark matter profile and a magnetic field _B M = 1 fj,G are adopted here. 



Under the assumption of equilibrium for the stellar component, we made use of the large wealth of available photo- 
metric studies to derive precise mass models for Draco. Under a general setup, we studied the preferred values for 
the dark matter halo length scale, its density normalization parameter as well as its anisotropy parameter. Results 
from numerical simulations of structure formation, together with a proper treatment of the effects of tides on the 
density profile, enabled us to correlate the concentration parameter and the initial virial mass of the dSph under 
consideration. In turn, this allowed us to further constrain the best-fit dark matter halo models for Draco. 

We then proceeded to an evaluation of the gamma-ray and electron/positron yield expected from Draco under the 
hypothesis that the dark matter is in the form of a pair-annihilating WIMP. To this extent, we resorted to illustrative 
cases of WIMPs of given mass and pair-annihilation cross section, dominantly annihilating into final states giving 
rise to the two extrema of a soft and a hard photon spectrum. For definiteness, and for illustrative purposes, we also 
considered theoretically well motivated benchmark supersymmetric models. 

We pointed out that unlike the case of the Milky Way galactic center, the spread in the estimate of the gamma-ray 
flux from Draco is significantly narrow, once the particle physics setup for the dark matter constituent is specified, 
and that Draco would appear as a point-like gamma-ray source in both ACTs and GLAST observations. In analogy 
with our procedure carried out in CPU2006 [lj for larger dark matter halos, we implemented a fully self-consistent 
propagation setup for positrons and electrons produced in WIMP pair annihilations, and we studied the subsequent 
generation of radiation in the radio frequencies from synchrotron emissions, and at higher frequencies from inverse 
Compton scattering off starlight and cosmic microwave background photons. 

We showed that, unlike in larger dark matter halos, as it is the case for the Coma cluster pj, in small, nearby 
objects the spatial diffusion of electrons and positrons plays a very significant role. As a consequence, the expected 
radio emission from Draco is spatially extended, and, depending upon the propagation setup and the values of the 
magnetic field in Draco, can provide a detectable signal for future radio telescopes. In some cases, we find that an 
extended radio emission could be detectable from Draco even if no gamma-ray source is identified by GLAST or by 
ACTs, making this technique the most promising search for dark matter signatures from the class of objects under 
consideration, i.e. nearby dwarf spheroidal galaxies. 

We finally showed that available data can accommodate the presence of a black hole in the center of Draco, even 
improving the fit to the data for some values of the black hole mass. The corresponding expected enhancement in 
the gamma-ray flux and in the radio surface brightness for cuspy dark matter halo profiles and an adiabatic growth 
of the black hole can be of several orders of magnitude. If the mass of the black hole is around or larger than 10 6 M©, 
WIMP models are predicted to give unmistakable astrophysical signatures both for future gamma-ray telescopes and 
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for future radio telescopes. 
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